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Abstract 

Let us denote by fi n the Birkhoff polytope ofnxn doubly-stochastic 
matrices. As the Birkhoff-von Neumann theorem famously states, 
the vertex set of fi n coincides with the set of all n x n permutation 
matrices. Here we seek a higher-dimensional analog of this basic fact. 

(2) 

Let fi n be the polytope which consists of all tristochastic arrays of 
order n. These are nxnx n arrays with nonnegative entries in which 

(2) 

every line sums to 1. What can be said about fi n s vertex set? 
It is well-known that an order-n Latin square may be viewed as a 
tristochastic array where every line contains n — 1 zeros and a single 

(2) 

1 entry. Indeed, every Latin square of order n is a vertex of fin , but 
as we show, such vertices constitute only a vanishingly small subset 

(2) 

of fi n s vertex set. More concretely, we show that the number of 
vertices of Clffl is at least {L n )^~°^ l \ where L n is the number of order 
n Latin squares. 

We also briefly consider similar problems concerning the polytope 
of n x n x n arrays where the entries in every coordinate hyperplane 
sum to 1. Several open questions are presented as well. 
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1 Introduction 



Let Q n C W 1 ' be the Birkhoff polytope, namely the set of order-n doubly 
stochastic matrices. The defining equations and inequalities of Q n are 

n 

Xij = 1 for all 1 < j < n 

t=l 
n 

Xij = 1 for all 1 < i < n 

3=1 

and 

> for all 1 < i,j < n. 
The vertex set of Q n is determined by the Birkhoff- von Neumann theorem 
[1, 7]. 

Theorem 1.1. The vertex set of Q n coincides with the set of permutation 
matrices of order n. 

We consider here some higher- dimensional analogs of the polytope Q n 
and ask about their vertex sets in light of Theorem 1.1. 

A line in an n x n x n array A is the set of entries obtained by fixing two 
indices and letting the third vary from 1 to n. A line of the form A(-,j, k) is 
called a column, a line of the form A(i, -, k) is a row and a line of the form 
A(i,j, ■) is a shaft. A coordinate hyperplane in A is the nxn matrix obtained 
by fixing one index and letting the other two vary. Such a hyperplane of the 
form A(-, ■, k) is called a layer of A. We denote the fc-th layer of A by A^. 
We denote the support of an array A by supp(A). 

(2) 

Let fin be the polytope of all tristochastic arrays of order n. Namely, 
nxnxn arrays with nonnegative entries in which every line sums to 1. Latin 
squares of order n can be viewed as two-dimensional permutations and it is 
easily verified that every Latin square of order n is a vertex of Vln^ ■ Does the 
natural analog of Theorem 1.1 hold true? As we show (Theorem 1.5), this is 
far from the truth. Of the v — v n vertices of Qn^ only fewer than t; 2 / 3 + ( 1 ) 
correspond to Latin squares. 

In section 3 we establish a similar phenomenon for a related polytope. 
Namely, now we consider nxnxn arrays of nonnegative reals in which 
every coordinate hyperplane sums to 1. Again it is shown that a natural, 
combinatorially defined set of vertices, are a vanishingly small subset of all 
vertices of this polytope. 
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1.1 Background material 



A Latin square L of order n is an n x n matrix with entries from [n] := 
{1, ...,n} such that each symbol appears exactly once in every row and col- 
umn. Equivalently, it is an n x n x n array A of zeros and ones in which every 
line has exactly one 1 entry. The correspondence between the two definitions 
is this: A(i,j, k) — 1 •<=> L(i,j) = k. We denote the number of order-n Latin 
squares by L n . 

The permanent of an n x n matrix A is defined as 



Per(A)=J2H< 



a-ei 



A lower bound on permanents of doubly stochastic matrices was conjec- 
tured by van der Waerden and proved by Falikman and by Egorychev [4, 3]. 

Theorem 1.2. If A is an n x n doubly stochastic matrix, then 

n) 

Per(A) > —. 



An upper bound on the permanent of zero/one matrices was conjectured 
by Mine and proved by Bregman [2]. 

Theorem 1.3. Let A be an n x n matrix of zeros and ones with ones in 
the i-th row. Then 

n 

Per {A) < Y[{ri\) 1/r '. 

i=l 

The following argument of van Lint and Wilson [6] utilizes these two 
bounds to derive an estimate for L n by constructing a Latin square A and 
bounding the number of ways to do this. Consider the n x n x n zero- 
one array representation of a Latin square layer by layer. Each layer is a 
permutation matrix, so that there are n\ choices for the first layer. Having 
already specified k — 1 layers, the number of choices for the /c-th layer can be 
expressed as the permanent of B, a zero/one matrix where bij = 1 iff a^t = 
for all k > t. Using the above upper and lower bounds on per(B) it follows 
that 
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Theorem 1.4. 

Ln= ((1+0(1))^) 



1.2 A higher dimensional Birkhoff polytope 
1.2.1 Definitions and motivating example 

Let Qn be the set of [n]^ 1 nonnegative arrays such that the sum of each 
line is 1. Thus, f2 n 1 ' ) = n , the set of order- n doubly stochastic matrices. 
Likewise, we call a member of fin a (d + l)-stochastic array. Maintaining 
the analogy, we let Sn be the set of [n] d+1 arrays of zeros and ones with a 
single one in each line. In other words, Sn consists of all (d + l)-stochastic 
arrays all of whose entries are zero or one. Thus, is the set of order n 

(2) 

permutation matrices and S n coincides with the set of order-n Latin squares. 
Members of Sn are called <i-permutations. 

We turn to investigate the vertex set of On • It is easily verified that 
every member of Sn is a vertex of fi n . However, as we show here Q$ 
can have numerous additional vertices. We find it instructive to present the 

(2) 

smallest such example. Namely, the following array A is a vertex of Q 3 . 
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To see that A is indeed a vertex, assume to the contrary that A = aB + 
(1 — a)C for some < a < 1 and B ^ C in Vtf \ If A(i,j, k) is or 1, then 
necessarily A(i,j,k) = B(i,j,k). So wherever A(i,j,k) ^ B(i,j,k), there 
holds A(i,j, k) = 2- 

Consider the graph G = G(A) whose vertices are the | entries of A, where 
two vertices are adjacent iff they are on the same line. Since A is tristochastic, 
it follows that B(i,j,k) + B(i',j',k') = 1 for every two neighbors (i,j,k) 
and (i',j',k') in G. Specifically, if B(i,j,k) = \ + e, then B(i',j',k') = 
\ — e. Consequently, the connected component of G which contains the 
vertices (i, j, k) and (i',f, k') is bipartite. The color of a vertex is determined 
according to whether the B entry is | ± e. However, it is easy to verify that 
G is connected and not bipartite, which proves our claim. 
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1.3 A scheme for constructing vertices 

(2) 

The above example suggests a construction for vertices of f2 n . Let A be an 
order-n tristochastic array whose support consists of exactly two | entries in 
each line. The graph G = G{A) defined as above is 3-regular and has 2n 2 

(2) 

vertices. As we now show, A is a vertex of Q n iff no connected component 
of G is bipartite. 

Indeed, suppose that G has a bipartite connected component with parts 
P and Q. Let A be the [n] 3 array with ±1 entries at the elements of P,Q 
respectively and everywhere else. Note that every line of A sums to zero. 
To see that A is not a vertex, note that A = X ~^ Y , where X,Y = A± |A are 
clearly tristochastic. 

Conversely, suppose that A = aB + (1 — a)C with 1 > a > and B ^ C 

(2) 

in f2 n is not a vertex. The same consideration that worked for the above 
example shows that the relevant component of G is bipartite. 

This discussion suggests that we construct A so that no connected com- 
ponent of G(A) is bipartite. This shouldn't be too hard, since G is 3-regular. 
Indeed, we suspect (but we still cannot show) that a typical tristochastic 
array with two |'s in each line is a vertex. This idea still yields the following 

(2) 

lower bound on the number of vertices of f2„ . 

(2) f— o(l) 

Theorem 1.5. The polytope f2„ has at least L% vertices. 

2 Proof of theorem 1.5 
2.1 The construction 

Let I = (ii,i2, ■ ■ ■ ,i n ) and J = ■ ■ ■ ,j n ) be two permutations of [n]. 

We let H(I, J) := {(z'i, Ji), (i 2 , ji), («2, J2), • • • , (i n ,jn), (kjn)} and call such 
a collection of index pairs an H-cycle. Note that H(I, J) = H(P, J') where 
I' = (h, h, ■ ■ ■ , in, h) and J' = (j 2 , js, ■ ■ ■ ,jn,ji)- Likewise, H(I, J) remains 
unchanged if we reverse the order of the i v and the j u . Consequently, the 
number of if-cycles is |n!(n — 1)!. 

Definition 2.1. Let n be an even integer. A double Latin square is an nxn 
matrix with entries from {1, |} where each symbol appears exactly twice 
in each row and column. 
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We say that a double Latin square X is Hamiltonian if the indices of the 
fc-entries of X constitute an i?-cycle for every fce{l,...,|}. (This explains 
the choice of the term if -cycle). 

Let A be a t x t matrix and let a 6 St be a permutation. We denote by 
(J (A) the matrix obtained from A by applying a to its rows. We need the 
following result from [5]: 

Proposition 2.2. Let A, B be two order | Latin squares and let a e §« be 

a cyclic permutation. Then the block matrix 

X -- 



A 


B 


a{A) 


B 



is a Hamiltonian double Latin square. 

It follows that the number of Hamiltonian order-??, double Latin squares 

is at least (f - 1)! • Ll = ((1 + o(l))^)^. 

We want to construct a tristochastic array A with exactly two |'s in each 
line, in such a way that G(A) is non-bipartite and connected (and therefore 
A is a vertex). 

The idea is to use a Hamiltonian double Latin square X to define the top 
| layers of A. We use the fact that X is Hamiltonian to complete A in such 
a way that G(A) is connected, and then "plant" an odd cycle in G(A) to 
ensure that G(A) isn't bipartite. 

Given a Hamiltonian double Latin square X, we use it as the "topograph- 
ical map" of the top | layers of A. Namely, A(i,j, k) = | X(i,j) = k. Let 
us observe the subgraph of G(A) spanned by the entries of A that reside in 
these top layers. Every positive entry A(i,j, k) — | comes from X(i,j) = k, 
and X has exactly two k entries in each line. Therefore this subgraph of 
G(A) is 2-regular. Moreover, since X is Hamiltonian, for every 1 < k < | 
the vertices of G(A) that correspond to supp(A fc ) constitute a cycle of length 
In. In other words, the subgraph of G(A) corresponding to the entries of the 
top half of A is the disjoint union of ~ cycles of length In. 

At this point, there are two | entries in every line that resides in one of 
the top | layers of A, and a single | entry in every shaft. 

We turn to define the next layer, As. +1 . Our purpose is to choose the | 
entries in this layer so as to form a single cycle of length 2n. The vertices 
of this subgraph should also be connected to each of the cycles in the top | 
layers. Clearly, if we manage to accomplish this task, then the part of G(A) 



(3 



that is already revealed is connected. Furthermore, note that every shaft 
contains a positive entry in the top half of A. Therefore, G(A) will remain 
connected regardless of our choices in the lower layers of A. 

In order to achieve our goals concerning A»+i, we want to find | index 
pairs (in, jn) such that X(i h ji) = I for all 1 < / < | and no two 

of them share a row or a column. We find such pairs successively as follows: 
Suppose that, for some k < |, we already have k pairs (h, ji), (ik, jk) 
with X(ii,ji) = 1, X(ik, jk) = k and no two pairs share a row or column. 
We claim that there is an additional pair (ik+i, jk+x) that does not share a 
row or column with any of the above index pairs, and X(ik+\, jk+i) = k + 1. 
Since X is a double Latin square, every row and column of X has exactly 
two elements that equal k + 1. Therefore at most 4k of these entries share 
a row or column with a previous pair. But 2n > 4k, so that such an index 
pair (ik+i, jk+i) must exist. 

We choose | more pairs of indices (in + ±, (in, in) in such a way 

that no two pairs of (ii,ji), (i n ,jn) share a row or a column. 

It is possible to rename, if necessary, the set of chosen pairs {(i a , i a )\a = 
1, . . . , n) as {(v, t v )\u = 1, . . . , n} for some permutation r £ £„. Let P be 
the permutation matrix of r. We next select a permutation a £ S n whose 
permutation matrix P' is such that P + P' consists of a single cycle. (We note 
that given r, there are exactly (n — 1)! possible choices for a). We achieve 
our aim by setting An +1 ■= \(P + P'). 

The purpose of our choices for A« +2 is to introduce an odd cycle into 
G(A). This odd cycle must use elements from the top half of A. Additionally, 
the indices of the | entries in v4™ +2 must avoid all index pairs used in As+i, 
so as not to create a shaft with three | entries. 

To this end, we seek two vertices x = (xi,x 2 , k) and y = (y\,y2, k) with 
x i 7^ Hi an d X2 7^ 2/2 that are connected by a path of odd length in the part 
of G(A) constructed so far. The construction of Ail + 2 will yield a length four 
path between x and y, ensuring that G(A) is not bipartite. This path will 
have the form x, x', w, y', y where x' = (xi,x 2 , f + 2), y' — (y±, y 2 , f + 2) and 
w is either (x±, y 2 , § + 2) or (y 1 , x 2 , | + 2). 

A simple counting argument shows the feasibility of this construction. 
Two vertices from the same layer can serve as x and y if their distance in 
that layer is odd and > 3. There are Q(n 2 ) such pairs in every layer with a 
total of Q(n 3 ) such candidate pairs. On the other hand, as we show below, 
only 0(n 2 ) such pairs are ruled out, so at least for large n a good choice of 
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such 2, y must exist. 

The reason that an entry cannot play the role of x is that its shaft meets 
supp(v4n +1 ). There are 0(ri) vertices in x's layer which might serve as y, and 
supp(y4n +1 ) has cardinality 2n, so only 0{n 2 ) pairs 2, y get ruled out for this 
reason. It remains to see how the pair x = (x\, X2, k) and y = (y±, y2, k) can 
be disqualified when both x's and y's shaft do not meet supp(An +1 ). This 
can happen only if both (21, y<i, f + 2) and (yx, X2, f + 2) are unavailable to 
us, namely A(xi, y%, | + 1) = A(yi, X2, § + 1) = \. There are only 0(n 2 ) such 
instances, one per each pair of vertices in the 2n-cycle residing in Ank+\. 

By doing these computations carefully, one shows that already for n > 10 
there must exist a good pair for the above argument. 

Next we need to complete supp(A™ +2 ). We are currently committed to 
three elements and 2n — 3 more ~ entries need to be chosen, so that altogether 
there are exactly two in each row and column. The locations that must not 
be chosen are those in the "shadow" of supp(A«. + i). It is easily seen that we 
need the following simple graph-theoretic claim. 

Proposition 2.3. Let G = (L, R, E) be a (n — 2) -regular bipartite graph with 
\R\ = \L\ = n > 6 and let M be a path of length 3 in G. Then there is a 
2-factor in G which contains the three edges of M. 

Proof. Let M = xi,X2,X3,x^. A bipartite graph with sides of size k and 
degrees > k/2 has a perfect matching. Let $ be a perfect matching in 
G \ {x\, X2, 23, £4}. Next let ^ be a perfect matching in G\{x2, 23} \ The 
desired 2-factor is $ U ^ U {(21,22), (22,23), (23,24)}. □ 

To recap, the graph G(A) is connected, it contains an odd cycle, and 
these properties are retained regardless of how the remaining § — 2 layers are 
completed. 

The remaining layers are constructed as follows. Let K be annxn matrix 
where K(i,j) = 1 or according to whether the shaft A(i,j,-) has one or 
two I entries. Each row and column of K has n — 4 one-entries. In other 
words, K is the adjacency matrix of an (n — 4)-regular bipartite graph which, 
therefore, has a 2-factor. This process can be completed layer by layer. This 
is just an existential argument and we next turn to estimate the number of 
ways in which our construction can be realized. 

To this end we will multiply the number of ways to construct the top 
half and the appropriate number for the bottom half. As stated above, there 

are L% + °^ ways to construct the top half. The estimate for the bottom 
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I — 2 layers is a slight variation on van Lint and Wilson's [6] approximate 
enumeration of Latin squares. By the van der Waerden bound [4, 3], a k- 
regular (n, n) bipartite graph H has at least ((1 + o(l)|) perfect matchings. 
By the same argument, there are at least ((1 + o(l)^-)" ways to complete a 
perfect matching in if to a 2-factor. The product of these two numbers is an 
overcount, since every cycle in the 2-factor can be split in two ways between 
the first and second 1-factors. Consequently, H has at least 

- + o(l)) Hk ~ 1 



e 2 V2 



2-factors. 

We think of K as the adjacency matrix of such an H, and each layer is just 
a 2-factor supported by K. With each choice, the edges of the chosen 2-factor 
are removed from H, which goes from being ci-regular to (d — 2)-regular. This 
yileds the following lower bound on the number of choices: 



n 



2<k<n-4, k is even 



o(l)) k{k 



The product of the bound for the top half and the bound for the bottom 
half yields a total of L„ °^' ) . 



3 A variation on the theme 

Here is another natural extension of the notion of doubly stochastic matrices. 
Namely, let to be the set of all arrays of nonnegative reals such 

that the entries in each coordinate hyperplane sum to one. The collection of 
such arrays clearly constitutes a convex polytope. Our goal is to investigate 
the vertex set of this polytope. Let us define Tn as the collection of all 
arrays of zeros and ones with a single one in each coordinate hyperplane. It 
is clear that Tn is included in the vertex set of £i . In view of the Birkhoff- 
von Neumann theorem it is natural to see how many of these vertices belong 
to Ti d) . 
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There is a natural bijection between tuples (oi, a^) G and members 
A G Tn which is given by <7i(i), ^(i)) — 1 for all 1 < i < n. In 
particular |T« | = 

As mentioned, every member of Ti^ is a vertex of £n , and we ask 
whether this polytope has any additional vertices. As it happens, such ver- 
tices are easy to construct. Here is the smallest example: 



At 



\ 



\ 

2 



(2) ~ 

Clearly A G We now consider the graph G(A) with vertex set 

supp(A) with an edge between every two vertices that lie in the same co- 
ordinate hyperplane. As in Section 1.2.1, we show that A is a vertex by 
observing that G(A) has no bipartite connected component. In the present 
case, G = K±. 

Our general construction is similar in nature to this example. We first 
construct an n x n matrix M with entries from [n] in which every row and 
column contains exactly two nonzero entries and where each integer in [n] 
appears exactly twice in M. We view M as a way to encode A as follows: 
M(i,j) = k for some k ^ says that A(i,j,k) = 1/2 and A(i,j,k') = 
for all k! ^ k. Also M(i,j) = means that A(i,j,l) = for all It is 
not hard to verify that if the graph G corresponding to M is connected and 
non-bipartite, then A is a vertex of T,nK 

We now turn to construct many such matrices M and thus generate many 
vertices for which are not in Tn \ Let 

H = (iijl), (*2, 32), ■ ■ ■ , (Wn). (hJn)} 

be an if -cycle and let 

M(ii,ji) = M(i 2 ,j 2 ) = 1 and M(i 2)]l ) = 2. 

The remaining entries of the if-cycle M(i a ,j a ) and M(i a+ i,j a ) are filled 
arbitrarily with the elements of the multiset {2, 3, 3, 4, 4, . . . , n, n}. Note 
that the resulting graph G is connected due to the if -cycle. It also contains 
the triangle 1), (12,32, 1)> (^2? ji, 2)}. There are |n!(n — 1)! choices for 

H and and ^ 2 2 "~^ ! ways to map the multiset to the nonzero entries of M. 

Altogether, this construction yields more than n! 4 > (t^j vertices of S^ 2 ' 1 . 
(2) 

It follows that Tn constitutes a vanishingly small subset of this vertex set. 
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4 Conjectures and some experimental results 



This paper raises many open questions. Here are several of them: 

(2) 

• Get a better estimate for the number of vertices of f2„ . 

• The analogous question for Qn with d > 2 seems completely open at 
this writing. 

(2) 

• The polytope f2„ is defined by requiring that one-dimensional subsets 

(2,) 

of the array sum to one. In the definition of En this is required of 
two-dimensional subsets. For larger d this suggests a whole range of 
possible polytopes to consider, depending on which sets of entries sum 
to 1. 

If we knew the support size of vertices in , we could make progress on 
these questions. By standard linear programming arguments, every vertex 
of fln^ has at least aff-dim(f^) zero coordinates. Since aff-dim(fii d) ) = 
(n— 1) , every vertex has support size at most n d+1 — (n— l) d+1 < (d+l)-n d . 

It follows that Q$ has at most Qj^-A^) < (i+r)^ 4 ^™ vertices. In 

particular, Qn^ has fewer than n 3 ™ 2 vertices. If we knew, say, that a typical 
vertex of Qn has support size < an 2 vertices, we could conclude that it has 
at most n ( 1 +°( 1 )) an vertices. 

We have conducted some numerical experiments to get a sense of the 
numbers. Using linear programming tools, it is possible to find the vertex 
that maximizes a randomly chosen linear objective function. Needless to 
say, this distribution on the vertices is by no means uniform. We neverthe- 
less hope that our experiments do tell us something meaningful about the 
properties of typical vertices. We selected the coordinates in the objective 
function independently from normal distribution. The average value of a in 
these experiments seems to increase slowly with n. We don't know whether 
the typical support size of a vertex converges to 3n 2 or to an 2 for some a < 3. 
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